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ABSTRACT 

This paper considers general relativistic (GR) effects in currently observed extrasolar 
planetary systems. Although GR corrections are small, they can compete with secular 
interactions in these systems and thereby play an important role. Specifically, some of 
the observed multiple planet systems are close to secular resonance, where the dynamics 
is extremely sensitive to GR corrections, and these systems can be used as laboratories 
to test general relativity. For the three-planet solar system Upsilon Andromedae, secular 
interaction theory implies an 80% probability of finding the system with its observed 
orbital elements if GR is correct, compared with only a 2% probability in the absence of 
GR. In the future, tighter constraints can be obtained with increased temporal coverage. 
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1. Introduction 

The perihelion advance of the planet Mercury is one of the classic tests of general relativity [1]. 
The recent discoveries of extrasolar planets [2] provide us with a new ensemble of solar systems to 
study periastron advance and other relativistic effects. The current sample of extrasolar planets 
[3] includes many orbits that are surprisingly close to the stars, with semimajor axes a « 0.05 AU 
and hence periods P ss 4 days. Such small semimajor axes imply that these planets experience 
larger GR effects than the planets in our solar system. One way to quantify the efficacy of GR is 
through the dimensionless parameter \x = GM t /(c 2 o). A planet in a 4 day orbit has a fi parameter 
nearly 10 times larger than that of Mercury, but direct measurement of periastron advance remains 
difficult because of the large number of observations required and because close planets tend to have 
nearly circular orbits. Fortunately, however, planetary systems with multiple planets and particular 
architectures allow for GR to exhibit much more pronounced effects. In these systems, secular 
interactions between the planets enforce time variations in the orbital elements (e.g., eccentricity 
e). These interactions are sensitive to the exact tuning of the system into secular resonance, and 
such tuning is affected by the relativistic corrections to the classical theory. As a result, these 
systems provide a new test of general relativity. In other systems, we find that this sensitivity to 
GR allows us to place new constraints on the system parameters. 

2. Theory of Secular Interactions 

We first outline the basic theory of secular interactions in multi-planet solar systems. In the absence 
of relativistic corrections, this topic has been widely discussed previously [4]. The basic result of 
these interactions is to force the orbits of participating planets to exchange angular momentum 
and thereby display time varying eccentricities (and other orbital elements). The amplitudes and 
time scales of these variations depend on the solar system configuration. Here we briefly outline 
the formalism and add the leading order relativistic correction. To the second order in eccentricity 
and inclination angle, the equations of motion for eccentricity ej and argument of periastron zuj 
decouple from those of inclination angle and the ascending node. Following standard convention 
[4] , we work in terms of the variables defined by 

hj = ej sin Wj and kj = ej cos Wj , (1) 

where the subscript refers to the jth planet in a solar system with TV" planets. The equations of 
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Vj = TljOj 



motion take the form 

dh i = J_dV 1 dkj_ = l_dVi^ (2) 

dt rija?j dkj dt njaj dhj ' 

where Vj the secular part of the disturbing function, rij is the mean motion, and aj is the semi- 
major axis (for the j'th planet). To consistent order in this approximation, the disturbing function 
can be written 

\ A n e ) + ^2 A jk eje k cos(wj - w k ) . (3) 

The physics of secular interactions is thus encapsulated in the N x N matrix A{j , where the matrix 
elements take the form 

rl ^ m k _ (i) Q GM*i 

A n = [4 L M^WWtt 01 *) + 3 ^-J ' (4) 

k=ij J ■> 

a 1 m k _ ( 2 ) . . . . 

A 3k = ~ n j 1 Mjf + — a jk a jkbl/ 2 (ttjfc) • (5) 

The quantities aj k are defined such that aj k = aj/a k (a k /aj) if aj < a k (a k < aj). The com- 
plementary quantities a,j k are defined so that a,j k = aj/a k = aj k if aj < a k , but a.j k = 1 for 
a k < aj. Finally, b^L(aj k ) and b^) 2 (aj k ) are the Laplace coefficients [4]. The diagonal matrix 
elements include the leading order correction for general relativity [5]. These corrections are small 
in an absolute sense, with /j~4x 10 -6 (aj/0.05AU) -1 , but they compete with the other terms and 
affect the eigenfrequencies when the system is near resonance. With the above definitions, the time 
variations of the eccentricity and argument of periastron are given by 

hj(t) = A ji sm(A^ + ft) , kj(t) = ^ Aji cos(A^ + ft) , (6) 

i i 

where the Aj are eigenvalues of the matrix A^j and the Aji are the corresponding eigenvectors. The 
phases ft and the normalization of the eigenvectors are determined by the initial conditions, i.e., 
the values of eccentricity ej and argument of periastron Wj for each planet at t = (taken to be 
the time when the orbital elements of the extrasolar planets are measured). 

3. General Relativity in Observed Systems 

To illustrate the action of GR in extrasolar planetary systems, we have calculated the secular in- 
teractions for two observed systems [3] with and without including the general relativistic terms. 
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We use the systems Upsilon Andromedae and HD160691 because they contain inner planets with 
a ~ 0.05 AU and outer planets with a ~ 1 AU. As shown below, this type of solar system architec- 
ture allows for GR to compete with secular interactions. Figure 1 shows the mean eccentricity (e) 
of the innermost planet, as driven by secular interactions and averaged over many cycles, plotted 
as a function of smi for the two systems. 

For Upsilon Andromedae (top panel), general relativity acts to damp the excitation of ec- 
centricity by secular interactions. Since the viewing angle is not measured, the resulting mean 
eccentricity (e) is shown as a function of sini. For large sini values, the inner planet would be 
driven to mean eccentricity values (e) ss 0.4 without GR, but only (e) « 0.016 when relativistic cor- 
rections are included. The observed eccentricity for the inner planet e f, s ~ 0.011 is much closer to 
the relativistic mean value. Since secular theory uses the observed orbital elements in the boundary 
conditions, the system always has some chance of displaying the observed (low) eccentricity of the 
inner planet (even if (e) ~ 0.4), so the implications for GR must be stated in terms of probabilities: 
If the observed eccentricity of the inner planet has a measurement error that is gaussian distributed 
with width a j, s = 0.015, then the probability of observing the system in its current state would 
be only 0.023 in the absence of GR. The probability of finding the inner planet with its observed 
eccentricity is 0.78 if GR is included. The validity of general relativity is thus strongly favored. 

For the HD160691 system (bottom panel), the inclusion of relativistic terms acts in the opposite 
direction, i.e., it leads to greater predicted values of (e). For values of sini ~ 1, the predicted mean 
driven eccentricities are small enough to be consistent with the observed low value e b s ~ 0. As sin % 
decreases, however, the level of eccentricity forcing increases and the system reaches a resonance 
near sini « 0.5. The observed low value of eccentricity, in conjunction with these considerations, 
thus constrain the viewing angle of the HD160691 system to be nearly edge-on; if we require 
(e) < 0.05, consistent with observational uncertainties, the viewing angle is confined to the range 
sini > 0.93 (i > 70°). Since inclination angles are notoriously difficult to measure in these systems, 
this constraint on i is quite valuable. For example, this limit on the viewing angle has important 
implications for the possibility of observing the inner planet in transit [6] . 

4. The Magnitude of Relativistic Effects 

Next we want to find an analytic criterion that characterizes the size of relativistic effects. In 
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these systems, general relativity is significant when the final term in equation (|4]) competes with 
the others. In practice, only the inner planet has significant relativistic corrections. We can also 
assume that the system is sufficiently hierarchical so that b^hfaji) = 30ji + 0{a 2 i), and that only 
one outer planet competes with the relativistic correction. In this limit, we obtain the requirement 



rrij 3 4GM* 

w^ 1 ~ 



where the subscript '1' denotes the inner planet and 'j' the outer planet. Both sides of this equation 
represent small dimensionless quantities. When their ratio is of order unity, however, relativistic 
effects compete with secular interactions, i.e., this constraint is equivalent to the requirement that 
one of the dimensionless fields II ~ 1, where 

= ±GMla) ^ , mj yl,M m y, aj y, at s-* 

c 2 m ja f ~ \mj) \M & ) VlAU/ V0.05AU/ ' K ) 

The second equality indicates that relativistic effects compete with secular interactions for a Jovian 
planet in a ~1 AU orbit perturbing an inner planet in a ~ 0.05 AU orbit. The relative size of the 
relativistic effect grows with increasing semi-major axis aj of the perturber, but the absolute size 
of the secular effect decreases. For closer perturbing planets (a,- <C 1 AU), interactions are stronger 
but relativity plays little role; for more distant planets {aj S> 1 AU), the interactions are weak but 
are dominated by relativity, which only makes the inner planet precess forward in its orbit. The 
condition for GR to compete with secular interactions can also be written in terms of time scales. 
GR itself defines a characteristic time scale r gr , 

" 2 " 3/2 ' a X \5/2/ M* \-3/2 



Tgr - 3(GM,)3/2 ~ 3011 y r U.05AuJ U.OmJ ' (9) 
the time required for the periastron to precess one radian forward in its orbit. When this time scale 
is comparable to the secular interaction time scale, r gr ~ T sec , then GR plays a significant role in 
the dynamics. For planets with ~ 4 day orbits, r gr is comparable to the secular interaction time 
scales for many of th.6 observed, extr&solcir planetary systems, where T sec — 10 2 - 10 5 yr [6]. 

5. The Sign of Relativistic Effects 

Figure 1 shows that GR can lead to either larger (HD160691) or smaller (Ups And) eccentricity 
forcing, compared with classical theory. However, most previous discussions of relativistic effects in 
planetary systems emphasize its stabilizing influence - the tendency for GR to reduce the amplitude 
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of eccentricity forcing. It is thus useful to explore when the two types of behavior occur, and what 
system parameters are required. 

For a two planet system, we consider an idealized case in which the inner planet has a nearly 
circular orbit, or at least cycles through the e = state. This condition is common in that close 
planets in observed systems tend to have nearly circular orbits [3]. In this limit, the eccentricity 
amplitude r/ of the inner planet (forced by the outer planet) can be written in the form 

V (Ai - A 2 )2 ' [W) 
where e 2 (o) is the initial eccentricity of the outer planet. For a two planet system, the difference in 
eigenvalues is positive definite and is given by the expression 

AA = Ai - A 2 = [(A u - A 22 f + 4A 12 A 21 } 1/2 . (11) 

Since the eigenvalues cannot be degenerate for a two planet system, secular interactions do not 
lead to resonance. Notice that only the diagonal matrix elements contain GR corrections and 
only the inner planet is close enough to the star for such corrections to matter. The first matrix 
element can thus be written An = Aq + 3ni/i, where Aq is the matrix element in the absence of 
GR and fi = GM*/(c 2 a>i) is the relativistic correction factor. With this construction, the question 
of whether GR acts to increase or decrease the amplitude of eccentricity oscillation depends on the 
sign of the derivative drj/dfj,, i.e., 



signf p-) = -sign(A + 3n lf i - A 22 ) = -sign 1 + H - ^(Zi) 1 /^ . ( 12 ) 



dfiJ 



mi /Oisi/2' 



m 2 a 2 



The second equality uses the approximations mi,m 2 <C M#, a\ <C a 2 , and ~ 3ai/a 2 ; the 
dimensionless field II is defined in equation ([8]) . The third term must dominate in order for the sign 
to be positive, and hence for relativity to lead to greater eccentricity amplitudes. Since a± < a 2 
by definition, relativity will lead to greater stability (smaller amplitudes) unless m\ > m 2 . As a 
result, relativity can amplify eccentricity forcing only for sufficiently massive inner planets, where 
the requirement for amplification can be written in the form 

mi > m 2 (^f 2 +4»M^) V2 . (13) 

This constraint thus implies that the inner planet mass m\ must be large compared to the outer 
planet mass m 2 . The two planet systems observed to date [3] show an interesting trend: None 
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of these solar systems have planetary masses that satisfy this inequality, and hence GR leads to 
greater stability in all of these (two planet) cases. 

Now we consider a three planet system in which the inner planet has little effect on the outer 
two planets due to its smaller mass and/or inner position. These characteristics apply to both 
HD160691 and UpsAnd. In this hierarchical limit, the 3x3 matrix Aij can be approximated by 
taking Ai 2 = A13 = A21 = A31 = 0. The eigenvalues of this reduced matrix then take the form 

Ai = A u , A 2 , 3 = \{{A 22 + ^33 ± [(A 33 - A 22 ) 2 + 4A 23 A 32 ] 1/2 } . (14) 

With this level of complexity, the system can have degenerate eigenvalues, e.g., when Ai = An is 
equal to either A2 or A3. Without loss of generality, suppose that Ai and A2 are nearly degenerate. 
The effect of GR is to add a small positive contribution to An and hence Ai (we again take An = 
Aq + 3nin), where the added term describes the relativistic precession of the inner planet's orbit. 
If Ai < A2, then GR brings the eigenvalues closer together and thus acts to increase eccentricity 
forcing. If Ai > A2, then GR makes the eigenvalues more unequal and acts to decrease eccentricity 
forcing. The full cubic equation for detL4 — XI] is more complicated, but allows for similar behavior. 

6. Conclusion 

For multiple planet systems with favorable architectures, the effects of general relativity can be 
significant (Fig. 1). When solar systems have two relatively massive outer planets and a third inner 
planet of smaller mass near secular resonance, GR effects are large enough to move the system in 
or out of a resonant condition. Since the resonance condition depends on planetary masses, which 
in turn depend on sini, and since small planets cannot survive in exact resonance, this effect can 
be used to constrain the allowed range of sinz in observed extrasolar planetary systems. For the 
HD160691 system, this constraint implies that sin i> 0.93. This GR effect thus provides another 
means of probing the properties of these solar systems. For other systems, general relativity causes 
the inner planet to precess forward in its orbit fast enough to compromise eccentricity excitation 
from a second planet, i.e., the mean eccentricity values are much smaller than they would be in flat 
space. This latter effect can be used as a new test of general relativity. For example, the Upsilon 
Andromedae system has an ~ 80% chance of being observed with its measured parameters if GR is 
correct, but only a ~ 2% chance if the GR corrections were absent. Although the weak-field limit 
of GR (including the perihelion advance of Mercury) is well established in our own solar system 
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[1] , the results of this paper add to our understanding by: [A] Developing a new relativistic effect 
- with amplified sensitivity - due to GR competing with secular interactions, and [B] Providing 
an independent confirmation of GR in other solar systems. In the future, additional extrasolar 
planetary systems and/or higher precision observational data can provide more stringent tests of 
general relativity. 
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Fig. 1. — The effects of general relativistic corrections on two extrasolar planetary systems. The mean 
eccentricity (e) of the innermost planet, as driven by secular interactions and averaged over many cycles, 
is plotted as a function of sinz for the Upsilon Andromedae system (top panel) and the HD 160691 system 
(bottom panel). Both systems have three planets detected to date. The predictions of secular theory are 
shown with relativistic corrections as the solid curves, and without relativistic corrections as the dashed 
curves. Notice that the relativistic terms act in opposite ways in the two systems: Inclusion of relativity acts 
to damp eccentricity excitation by the secular interactions in the Upsilon Andromedae system. In HD160691, 
however, they allow the system to approach a resonance for sini « 0.5. 



